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Multi-user detection method 



The present invention concerns in general terms a multi-user detection method 
(MUD). More precisely, the present invention relates to an iterative method of 
eliminating interference between users (Multiple Access Interference or MAI). The 
present invention applies more particularly to mobile telephony in DS-CDMA mode 
5 (Direct Sequence - Code Division Multiple Access) that is to say to mobile telephony 
using a code division access mode with spectral spreading by direct sequences. 

In a DS-CDMA mobile telephony system, the different users are separated by 
multiplying each symbol of the user by a spreading sequence peculiar to it, also 
referred to for this reason as the signature of the user, the different signatures ideally 

10 being chosen so as to be orthogonal. The spreading sequence frequency (chip rate) 
being greater than the frequency of the symbols, the signal transmitted by each user is 
distributed (or spread) in the frequency space. On reception, the signal of a user is 
separated by virtue of a filtering adapted to the corresponding signature. This filtering 
is also referred to as "despreading". The ratio between the band occupied by the 

15 spread signal and the band occupied by the information signal is referred to as the 
spread factor. 

The signatures employed must have good correlation properties, namely a very 
pronounced auto-correlation peak and low intercorrelation values. 

The first of these two characteristics allows synchronisation of the sequence 

20 received. It is very useful when the transmission channel of a user has several 
propagation paths. This is because each path must then be isolated by virtue of a 
filtering adapted to the signature and to the delay of the path. It is possible to take 
advantage of the diversity of propagation within the channel to increase the signal to 
noise ratio on reception. To do this, a bank of adapted filters, separating the different 

25 paths, is used and the outputs of them are combined. The most widespread 
combination is MRC (Maximum Ratio Combining), which consists of multiplying the 
signal output from each adapted filter by the conjugate of the complex multiplicative 
coefficient introduced by the channel on the path concerned. The resulting filtering 
operation is a filtering adapted to the equivalent filter of the channel. Through its 

30 structure, the receiver thus formed is referred to as a rake receiver. Naturally, perfect 
separation of the paths takes place only if the auto-correlation is a Dirac. In practice, 
however, separation is not complete and leaves multipath interference which is also 
referred to as self-noise. Fig. 1 depicts schematically a DS-CDMA system with K 
users. The data of a user k are spread in frequency by the corresponding signature in 
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the module 100 k before being transmitted over the channel 110k having P paths. On 
reception, for a given user k, the signals being propagated according to the different 
paths p=L.P of the channel are separated by adapted filters 120 ki i...l20k,p (only the 
battery of filters of the user k has been shown) before being weighted by a set of 
5 complex coefficients ck p . The signals thus weighted are summed (140k) and the 
resulting sum output from the rake receiver is subsequently detected to provide an 
estimation of the data of the user L In the case of a downlink (links from a base 
station to a mobile terminal) the channels 1 to K are identical, whilst they are 
different in the uplink (links from mobile terminals to the base station). The first case 

1 0 can, from this point of view, be considered to be a particular case of the second. 

The second characteristic set out above guarantees a low level of interference 
between two distinct users. Nevertheless, here also, in practice, the intercorrelation 
between two signatures is rarely zero. This is notably the case in a so-called dazzle 
situation (the Near-far effect) where a high-power signal received from a user 

15 interferes with the reception of a low-power signal coming from another user. 
Moreover, when the number of users is high, close to the spread factor, the sum of the 
interferences of the different users, low if taken in isolation, can have very disturbing 
effects on detection. 

In order to combat the multi-user interference, several methods have been 
20 proposed. A review thereof will be found in the article by Simon Moshavi entitled 
"Multi-user detection for DS-CDMA communications", which appeared in IEEE 
Communications Magazine, Oct. 1996, pages 124-136. Amongst existing multi-user 
techniques, the techniques of subtractive elimination (Subtractive Interference 
Cancellation) have good performance with reasonable complexity in use. The general 
25 idea of it is simple; from a first detection at the output of the adapted filter, the 
contributions to the interference suffered by the other users are reconstructed by 
respreading. This interference is then subtracted from the signal received in order to 
provide a cleaned signal at a subsequent detection step. According to the way in 
which the subtraction is carried out, there may be parallel elimination (PIC, standing 
30 for Parallel Interference Cancellation) and serial elimination of interference (SIC, 
standing for Serial Interference Cancellation). 

Let it be assumed first of all that the signal of each user propagates on a single 
path as far as the receiver. 
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The parallel elimination method is illustrated in Fig. 2: the received signal is 
filtered by a battery of adapted filters (200i,200 2 ,..,200k), each adapted filter relating 
to a given user. After detection (210 k ), the estimated symbols are respread (220 k ) 
spectrally by means of the signature of the user in question before being filtered by a 
5 filter modelling the transmission channel (230*). There is thus, at the output of (230*), 
an estimation on the contributory share of the signal received attributable to the user L 
The sum of the contributory parts of the other users is then subtracted (at (240*)) in 
order to obtain a cleaned signal S* (I \ This cleaned signal can directly be the subject of 
detection after despreading or the elimination process can be reiterated. The detection 

10 being of better quality at each iteration, there is then obtained, at the end of successive 
iterations, signals S* w which are more and more rid of the multi-user interference. 

The serial elimination method is illustrated in Fig. 3: the signals received by the 
different users are first of all ordered in decreasing order of power, that is to say l 9 . my K. 
The procedure then consists of successive eliminations of the contributory shares, 

15 commencing with the signal of highest power. To this end, the SIC detector has a 
series of stages in cascade, each eliminating the interference due to a particular user. 
The first stage works on the antenna signal and each subsequent stage receives as an 
input the output of the previous stage. Each stage has an adapted filter (300*), a 
detector (310 k ), a module (320*) for respreading the symbols, a filter (330*) modelling 

20 the transmission channel k and a subtractor (340*) eliminating the contribution due to 
the user k. Each stage also supplies at the output of the detector (310*) a decision on 
the symbol received, S*, and the interference elimination process ends at the K th stage. 

The techniques set out above apply well to the simple situation where the 
transmission channel of a user has a single path. In this case, the filter modelling the 

25 channel can be limited to multiplication by a complex coefficient. When the channels 
are multipath, the situation is on the other hand much more complex since it is 
necessary to eliminate both the multipath interference and the multi-user interference. 
An iterative detector with subtractive elimination of the multi-user interference when 
there are multiple paths was proposed in an article by M.C. Reed et al. entitled 

30 "Iterative Multiuser detection using antenna arrays and FEC on multipath channels" 
published in the IEEE Journal on Selected Areas in Communications, Vol. 17, N° 12, 
Dec. 1999, pages 2082-2089. Each iteration of the detection comprises an adapted 
filtering, a channel formation and a combination of the rake type. The method 
proposed presupposes however that the attenuation coefficients, the phase rotations 
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and the directions of anival of all the paths of all the users are determined. This 
determination is carried out externally to the detector,, prior to the sequence off 
iterations and therefore on signals interfered with by multi-user and multipath 
• interference. As a result the elimination of this interference is necessarily 
5 approximate. 

A first aim of the invention is to propose a detection method and device which 
allow elimination of multi-user interference when there are multipath channels 
without requiring an explicit calculation of the attenuation coefficients, phase 
rotations and directions of arrival of the paths of the different users. 

10 A second aim of the present invention is to propose a detection method and 

device which improve the elimination of multi-user and multipath interference. 

The object of the invention is defined by a multi-user detection method with 
elimination of interference between users, each user transmitting modulated data in 
the form of symbols on a transmission channel, each transmission channel 

15 comprising at least one propagation path and each propagation path arriving at an 
array of reception antennae according to a direction of arrival, the method 
comprising at least one sequence of steps for each user, each sequence comprising: 

(a) a reception step decomposing each antenna signal into filtered signals 
issuing from the different paths of the said user and combining the said filtered 

20 signals by means of a first plurality of coefficients in order to form an estimation of 
the signal transmitted by the user; 

(b) a step of estimating the contribution of the user to the signals received by 
the different antennae from the said estimation of the signal transmitted and a second 
plurality of coefficients obtained from the said first plurality of coefficients; 

25 (c) a step of eliminating interference from the cleaned antennae signals; 

the cleaned antennae signals supplied by at least one first sequence being used 
as antennae signals by at least one second sequence. 

According to a first embodiment, the first plurality of coefficients comprises a 
first set of complex coefficients and a second set of complex coefficients and the 
30 filtered signals are subjected to a channel formation step in order to form signals of 
paths by means of the said first set, the said path signals then being linearly 
combined by means of the said second set in order to supply the said estimation of 
the signal transmitted, the coefficients of the first set being adapted so as to minimise 
a plurality of first error signals between a reference value of the transmitted signal 
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and the said path signals, the coefficients of the said second set being adapted so as 
to minimise a second error signal between the said estimation and the said reference 
value. 

According to one characteristic of the invention, the said second plurality of 
5 coefficients comprises a first set of complex coefficients and a second set of complex 
coefficients, the coefficients in the said first set of the second plurality being obtained 
from the arguments of the coefficients of the first set of the first plurality and the 
coefficients of the said second set of the second plurality being obtained from 
coefficients of the said second set of the first plurality. 
10 Advantageously, the coefficients in the said first set of the second plurality are 

obtained from a linear regression on the arguments of the coefficients of the first set of 
the first plurality. 

According to a second embodiment of the invention, the first plurality of 
coefficients consists of a set of complex coefficients and the filtered signals are 
15 linearly combined by means of the said set in order to supply the said estimation of 
the signal transmitted, the coefficients of the said set being adapted so as to minimise 
an error signal between the said estimation and a reference value. 

According to one characteristic of the invention, the said second plurality 
( wc,i,k ,wu) of coefficients comprises a first set of complex coefficients (hv,*) and a 
20 second set of complex coefficients O,*), the coefficients (m*) of the said second set 
of the second plurality being obtained by: wjk=gi*lgk where gi,k is an estimation of 

the norm of the sub-vector ba= . , L being the number of antennae in the array, 

where gk is a mean of the g& values on the different paths, the coefficients ( vw,,-,* ) of 
the said first set of the second plurality being obtained from the arguments of the 

25 coefficients ( bejjt ) of the said set of the first plurality. 

Advantageously, each sequence comprises a step of estimating the symbols 
transmitted from the estimation of the signal transmitted in order to obtain first 
estimated signals, a step of demodulating the first estimated symbols in order to 
obtain estimated data, a step of channel decoding of the said estimated data followed 

30 by a channel recoding and a remodulation in order to obtain second estimated 
symbols. 
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According to a second variant of the invention, the users are classified in order 
of received power and the interference is eliminated by subtracting one after another 
the contributions of the different users, commencing with the users with the highest 
received powers. 

5 Advantageously, the estimations of the transmitted signals of the users being 

considered to be an estimations vector with K components where K is the number of 
users, the said vector is subjected to a transverse matrix filtering. 

In addition, the estimated symbols of the users being considered to be a 
symbols vector with K components, the said symbols vector can be subjected to a 
10 postcursor matrix filtering and the output of this filtering is subtracted, vector by 
vector, from the output of the transverse matrix filtering. 

The invention is also defined by a multi-user detection device comprising 
means adapted to implement the method disclosed above. 

The characteristics of the invention mentioned above, as well as others, will 
15 emerge more clearly from a reading of the description given in relation to the 
accompanying figures, amongst which: 

Fig. 1 depicts schematically a DS-CDMA telecommunications system; 
Fig. 2 depicts schematically a multi-user detection device of the PIC type; 
20 Fig. 3 depicts schematically a multi-user detection device of the SIC type; 

Fig. 4 depicts a first receiver of the rake type useful to the device according to 
the invention; 

Fig. 5 depicts a second receiver of the rake type useful to the device according 
to the invention; 

25 Fig. 6 depicts a multi-user detection device according to a first embodiment of 

the invention; 

Fig. 7 depicts a multi-user detection device according to a second embodiment 
of the invention; 

Fig. 8 a depicts a multi-user detection device according to a third embodiment of 
30 the invention; 

Fig. 8b depicts in an equivalent fashion a multi-user detection device according 
to the third embodiment of the invention; 

Fig. 9 depicts a multi-user detection device according to a fourth embodiment of 
the invention. 
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The general idea at the basis of the invention is to propose an adaptive reception 
method of the rake type, fitted into a multi-user detection process. The interference 
rejection process benefits from increased precision by virtue of the estimation of the 
5 respective contributions of the different users in the received signal. Conversely, the 
process of adapting the complex coefficients used in the reception method of the rake 
type benefits from the rejection of the multi-user interference by operating on the 
cleaned signals. 

10 According to a first variant, the invention uses a receiver of the rake type as 

illustrated in Fig. 4. 

The receiver depicted is dedicated to the reception of symbols transmitted by a 
user k. The system has L antennae represented schematically by the block (400). 
Each of the L antenna signals is then correlated by a battery of P adapted filters, where 
15 P is the number of propagation paths of the channel, each filter being adapted to a 
given path i. The set of L.P adapted filters is represented diagrammatically by the 
block (41 0 k ). 

If the notation is given to the contribution received by the antenna £ y 
£=l,...,L attributable to a path /, i=l„P of the channel k, k=l..K it is possible to 
20 write: 

jc^(0^u(/).e y ^^*(^H"«.w(0 (1) 

with 

25 £^(0w(0H^-i)-M^ ( 2 ) 

where Sk(t) is the symbol sent by the user k, nu*(f) is the assumed white 
additive Gaussian noise at the antenna £ for the path i of the user k y including the 
thermal noise, the multipath interference and the multi-user interference, d is the pitch 
of the array and X is the wavelength. 

30 The propagation along the path i of a signal sent by the user k can thus be 

modelled by an attenuation coefficient au(0 a phase rotation vu(/) . The total 
phase shift £w.*(/) results on the one hand from the phase rotation vu(t) induced on 
the path i and on the other hand the difference in operation between the antenna £ and 
the reference antenna 1 at the angle of arrival &.*(*) . 
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The L.P adapted filters supply the signals xtj* , £=l.Z, i=l..P, the signals 
issuing from the filters adapted to the same path i being directed to a channel former 
associated with this path. The P channel formers are depicted in the form of a single 
block (420k). The output signals from the channel formers, y,\k , are written: 

yuk =^Txejjt. bej t k 



(3) 



where bujt are complex weighting coefficients. Operation (3) can also be 
written in an equivalent manner in the following vector form: 

( xxjjt(t/\ 



X2,i,k(t) 
\XL.ik(t)) 



(3') 
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where bu.k(t)=ae^k,cxp(-j y/tu) with £ = (i,.„,L) 

Secondly, the output signals yuk , that is to say the different paths of the signal 
from the user k, are combined by the module (440k) in order to supply a combination 
signal zt . The combination signal z* can then be written: 



zk(t)=Ck(t) T .yk(t)4ajc(t)ci,k(t) ■ cp,k(tj). 



(4) 



where ca(%)=p/,*.exp(— £u) with i = (l,...,P) are complex weighting coefficients. 

The coefficients be,i.k are adapted so as to minimise, in root mean square value, 
the error signals expressed by: 

dukOhqUthy^t) (5) 
Likewise, the coefficients a,k are adapted so as to minimise, in root mean square 

t 

value, the error signals expressed by: 

t?k(t)=q\(t)-zk(t) (6) 
where q\k(t) and q n k(t) are reference values corresponding to the observations output 
from (420 k ) and (440 k ) respectively. The reference values q\k(t) and q"k(i) are 
chosen so as to be equal: 

q\k(t)=q n k{t)=qk(t) where q^t) is a common reference value which does not depend on 
the path /. 

In order to adapt the coefficients and ajk , it is possible to use, as a common 
reference value qk(t), either directly pilot symbols Pk(t) 9 or estimated symbols, before 
or after channel decoding, as will be seen below. Two possibilities of reference values 
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can also follow each other, the adaptation of the coefficients bejjt and a,* then starting 
with pilot symbols and continuing with estimated symbols. 

Whatever the case, during the adaptation, the coefficients are modified so as to 
minimise the mean values E^kf) and E^'kf) . Because of the feedback of the error 
5 signal f'u, the values y/tu converge progressively towards the values 

and consequently compensate for the phase differences due to the channel and to the 
differences in operation between antennae. As a result the channel formers are not 
content to form beams in the direction of the different paths but also at least partially 
compensate for the phase rotations caused by the different paths. Such a type of 

10 channel former will be referred to hereinafter as a "channel former with phase 
precompensation". The module (440 k ) then effects a combination on the signals of 
paths precompensated in phase issuing from the channel formers and corrects the 
phase shift remainders. 

The coefficient adaptation algorithm is based on a minimisation, in root mean 

15 square value, of an error signal. Many algorithms have been proposed for this purpose 
in the literature and a description of the main ones can be found in the journal article 
by L.C. Godara entitled "Application of antenna arrays to mobile communications, 
Part II: Beam-forming and direction-of-arrival considerations", published in 
Proceedings of the IEEE, Vol. 85, N° 8, August 1997, pages 1 195-1245. 

20 For example, it is possible to use an algorithm of the LMS (Least Mean Square) 

type or an algorithm of the RLS (Recursive Least Square) type as disclosed in the 
French application entitled "Bi-modular adaptive CDMA receiver" filed in the name 
of the applicant. 

The receiver described above will be referred to hereinafter as a "receiver of the 
25 first type". 

According to a second variant, the invention uses a receiver of the rake type as 
illustrated in Fig. 5. 

The receiver depicted is dedicated to the reception of symbols transmitted over a 
30 channel k. The system includes L antennae represented schematically by the block 
(500). Each of the L antenna signals is then correlated by a battery of P adapted filters 
where P is the number of propagation paths of the channel, each filter being adapted 
to a path L The set of L.P adapted filters is represented schematically by the block 
(510 k ), the L.P output signals, xw, £=1.JL 9 i=l..P of the adapted filters being 
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transmitted to a module linearly combining these output signals by means of complex 
coefficients ht t u represented by the block (520k) in order to supply a combination 
signal zk . The combination module (520k) performs the following operations: 

P L 

Zk = /^ X£,i,k~bej,k (7) 

where the values beM(t)=aejt.exp(~j y/tu) with £=l,...,L, i=l..P, are complex 
weighting coefficients. Equation (4) is written in an equivalent manner: 



Zk(t)=bk(t) T 'Xk(t)Mh,\,k(t)bixk(t} • J>Lxk(t)..h,p,k(t)bi,p,k(t} • i>L t p,k(tj). 



X2Xk(t) 
XL,\ t k(t) 

xi,p*(i) 



(V) 



^xlM*)) 

The coefficients bejjt are adapted so as to minimise, in root mean square value, 
the error signal expressed by: 
10 £k(t)=qk(t)-zk(t) (8) 

where qk(t) is a reference value corresponding to the observed signal zk(t) . 
In order to adapt the coefficients it is possible to use, as the reference value, 
either directly pilot symbols p k (t), or estimated symbols, before or after channel 
decoding, as will be seen below. The two possibilities of reference values can also 
15 follow each other, the adaptation of the coefficients &,/,* then starting with the pilot 
symbols and continuing with estimated symbols. 

During the adaptation, the coefficients are modified so as to minimise the 
mean value i?(|a| 2 ) . Because of the feedback of the error signal at , the values y/m 
converge progressively towards the values and consequently 

20 compensate for the phase differences due to the channel and to the differences in 
operation between antennae. As a result the module (520k) effects both a channel 
formation and a compensation for the phase rotations caused by the different paths. In 
addition, the module (520 k ) effects a real weighting of the contributions of the 
different paths of each user. 
25 The adaptation algorithm can, there too, be an algorithm of the LMS (Least 

Mean Square) type or an algorithm of the RLS (Recursive Least Square) type as 
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disclosed in the French application entitled "Single-module adaptive CDMA receiver" 
filed in the name of the applicant. 

The receiver described above will be designed hereinafter as a "receiver of the 
second type". 

5 

Fig. 6 depicts a multi-user detection device according to a first example 
embodiment of the invention. The multi-user interference is eliminated therein 
according to a parallel mode. 

The outputs of the antennae I..L in the array are all connected to a battery of 

10 receivers (600i,.. ,600k) of the first type or second type, as illustrated respectively in 
Figs. 4 and 5. The signals obtained at the output of these receivers are then subjected 
to a sequence of processing steps. This sequence comprises first of all a decision step 
(610i,..,610k) supplying in soft or hard form an estimation of the symbols transmitted 
by the different users, that is to say Sk(t) . If the decisions on the symbols taken at 

15 (610k) have soft values, the estimated symbols Sk(t) are expressed as a weighting of 
the symbols of the constellation by their respective a posteriori probabilities. If the 
decisions on the symbols taken at (610 k ) have hard values, the estimated symbols Sk(t) 
result from a thresholding of the observed signal z k (t) followed, possibly, by a 
weighting by a factor determined a priori. The estimated symbols are then 

20 demodulated (620i,. ..,620k) in order to supply estimated data dk(t). The estimated 
data are subjected to a step of deinterleaving and channel decoding (630i,. . .,630 K ) and 
then recoded and interleaved once again (640i,...,640k) in order to supply a more 
reliable estimation of the coded data. Advantageously, the decoding operation is 
performed by means of a decoding with soft inputs and outputs (SISO, standing for 

25 Soft Input Soft Output). It should be noted that recoding is not necessary if the 
encoded data are recovered directly at the decoding output. For example, if use is 
made at (630i,...,630 K ) of an algorithm of the SOVA (Soft Output Viterbi Algorithm) 
type the coded data can be recovered directly along the maximum likelihood path. 
The newly estimated data ctk{i) are remodulated (65 0 k ) and the symbols obtained 

30 s\{i) , are respread spectrally (660 k ) by the signal of the user L Naturally, if the 
decoding used is of the soft value type, remodulation is also "soft" in the sense that it 
supplies symbols resulting from the weighting of symbols of the constellation, 
weighted by their respective a posteriori probabilities. Each respread symbol s\{t) is 
then filtered by filter (670 k ) whose function will be explained below, whose role is to 
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model the effect of the transmission channel k and the differences in operation 
between antennae. There is thus, at the output of the filter (670k), an estimation of the 
"contribution" of the user k to the L antenna signals, the said contribution being in the 
form of an L-tuplet of signals (xe,k)t=i.x • 
5 There is deducted at (680*), for each user k, from the L-tuplet of signals (xt) 

received by the L antennae, the L-tuplet (£ Xe #) in order to eliminate the 

contribution due to the other users. In this way K L-tuplets of antenna signals are 
obtained, with the multi-user interference partly removed. The elimination process 
can be iterated by repeating the sequence (600k)..(670k) as from these signals and 

10 deducting, at each end of an iteration, from the signals (xt) received by the antennae, 
the newly calculated contributions. Each sequence of a user of the second iteration 
and the following iterations operates on a set (L-tuplet) of antennae where the 
contributions of the other users were eliminated at the previous iteration. The 
estimation of the symbols transmitted by the different users thus gains in reliability as 

1 5 the iterations progress. 

It should be noted that the demodulation, deinterleaving, channel 
decoding/coding, interleaving and remodulation operations are optional (the blocks 
620k, 630k, 640k, 650k are then absent). The spectral respreading can be applied 
directly to the estimated symbol $k(t) . However, the use of a channel decoding (with 

20 of course corresponding channel coding on sending) reduces the error rate and 
therefore improves the reliability of the estimation of the interference and, 
consequently, the reliability of the transmitted symbols. 

The structure of the filter (670 k ) depends on the type of receiver used at (600 k ). 
25 The two types of receiver envisaged will be examined successively: 

If the filter (600 k ) is of the first type, the filter (670 k ) performs the following 
operation: 

30 jca(/)=]f]vw,u(/).w40i^-^) ( 9 ) 
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where the fa values are the estimated delay times of the different paths i of the 
channel k obtained from the time differences between the correlation peaks calculated 
at (410k); 

wjt are complex weighting coefficients modelling the "apparent" channel, that is 
5 to say where the different paths have been precompensated in phase, seen by the path 
recombination module (440k); 

mj* are complex weighting coefficients modelling the differences in operation 
between antennae and the phase precompensation of the module (420 k ). 

10 According to a first variant embodiment of the invention, the coefficients wjk 

are chosen such that wa=c? ft , in a similar manner to the MRC algorithm applied to the 

path recombination. Thus multiplicative terms of the channel corresponding to the 
attenuations and the phase difference remainders applied in (440 k ) are reconstructed. 
According to a second embodiment, it is considered that the phase difference 
15 remainders are negligible and therefore that the coefficients a* are essentially real. 
There is then chosen: 

Mu=jcu| or ugt=cu or wu=Re[ck\*] (10) 
where Re[J represents the real part. 
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In the same way, the coefficients bujt are considered to be pure phase 
differences. The terms mjjc can then be chosen so as to reconstruct the phase 
differences corresponding to the coefficients but applied in (420 k ). 

According to a first variant embodiment, there is then chosen: 
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According to a second variant of the invention, a linear regression is effected on 
the values y/tj* as described in French application FR-A-0011160 filed on 29.8.2000 
in the name of the applicant and incorporated here by reference. More precisely, for a 
30 given path z, the values n* and (jx* which minimise the mean square deviation J are 
sought: 

J( vw,<p,\k)=^(vi t k +(£-1 ). <pij -y/ej,k J (12) 

These two values are obtained by cancellation of the partial derivatives of the function 
J , and there is then obtained: 
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vu ~if/i t k ~(pi,k. (£-1) 



(13) 



with ^ijt^^jftM and £ - 1 = j £ {£ - 1) 
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and fa Jt-V-W^v *£l 

(£-iy-£-i 



• (14) 



where X designates the average of Xover the L antennae. By means of equation 
(13), it is then possible to calculate v,\* . From the values n* and <pi* there is then 
chosen: 



w#,i\*=exp(/( vu+(^-l ). <pu)) 



(15) 



which gives a less noisy estimation than the one supplied by equation (1 1). 

15 If now the filter (600 k ) is of the second type, the filter (670 k ) performs the 

following operations: 



p _ 



(16) 



20 The coefficients uu are obtained by considering that the energy of the radiation 

diagram of the channel former relating to a path is proportional to the attenuation 
undergone by the signal along this path. If there is noted: 
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, ^A==||bu|| and gk 



there is then chosen uuk =gu /gk ( 1 7) 

The norm can, for example, be approximated by £a=maxj£a*j or by 

- ^ 

gt\k —biu =1/ L . 
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The coefficients vw.a are obtained, for their part, as before by equation (1 1): 
VMf 1 a=exp(/^,/ 1 *) where ^w;*=Aictan(lii][&w;jt3/Re[Z?^i\jt]) with possibly elimination of 
the phase noise by linear regression, as indicated by equations (12) to (15). 

5 When a filter of the first or second type is used, the coefficients but and, where 

applicable, ca , must be initialised at each iteration. 

For the first iteration, the coefficients will be initialised by the values: 

for a receiver of the first type: 
10 buk(0)=S(£-£o) 9 Vi and o*(0)=e,Vi (18) 

where S is the Dirac symbol, £o is an antenna number and c a given complex 
coefficient; 

15 for a receiver of the second type: 

bu*(0}=bJ(£-to) 9 Vi (19) 

where S is the Dirac symbol, £a is an antenna number and b a given complex 
20 coefficient. 

In both cases, the initial configuration chosen is a multipath omni-directional 
configuration. 

Alternatively, if there is an initial estimation of the attenuation coefficients, 
25 directions of arrival and phase rotations of the different paths, it is possible to initialise 
the coefficients, at the first iteration, by the values: 

for a receiver of the first type: 

30 W(0)-expH(va(0)+2^M.cos^(0).(^-l))) and a*(0)=&*(0) (20) 

for a receiver of the second type: 



fc,*(0)^a(0).exp(-yX ^ 



(21) 
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where ^u(O) , v;,*(0) , &jt(0) are respectively the initial estimations of the 
directions of arrival, phase rotations and coefficients of attenuation for the different 
paths. 

5 At the second iteration and at the subsequent iterations the coefficients but and, 

where applicable, a* are initialised by the values of the coefficients obtained at the 
end of adaptation at the previous iteration. More precisely, if it is agreed to initialise 
the time to zero at the start of each new iteration and if the length of the block of 
symbols processed by an iteration is denoted T, there is chosen: 

10 

b {n \M0)=b (n - ]) eM(T) (22) 



where the index n indicates the number of the iteration. 

If the adaptation process has not converged at the previous iteration, that is to 
15 say, for example, if the modules of the error signal or signals are greater than a 
predetermined threshold during this iteration, it may be preferred to reinitialise the 
said coefficients to values given by (18), (19) or (20), (21). It will also be possible to 
use as initial values a mean, possibly weighted by a forgetting factor S , of the final 
valuations of the previous iterations, that is to say: 

20 

b {n \jm=^ m b {n - m \jAT) (23) 

m=\ 



After the phase of initialising the coefficients, it is important to specify what 
reference values qk(i) are used by the receivers (600 k ): it is possible to use pilot 
25 symbols or symbols Sk obtained after soft or hard decision at the output of (610 k ) or 
symbols obtained after channel decoding/recoding at the output of (650k). In order not 
to delay the feedback of the error signal by waiting for operations (610k) to (650 k ), the 
values of the symbols obtained at the previous iteration will advantageously be used, 
that is to say: 

30 

q {n \(t)=s«"- l \(t) (24) 



where q (n \(t) designates the reference value used at time / of iteration n. 
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It can also be envisaged taking as the reference value a weighted mean between 
a symbol output from (610 k ) of the current iteration and a corresponding symbol 
output from (650k) of the previous iteration, that is to say: 

5 ^ <B) A(0= 7 7^ ,(/, ~ 1) j fc (O+Cl-^-^^CO where tj is a weighting coefficient ( 0 < 77 < 1 ). 

Fig, 7 illustrates a multi-user detection device according to a second example 
embodiment of the invention. 

Unlike Fig. 6, the elimination of the interference is performed here in a series of 

10 K steps, K being the number of users. The contributions of the different users to the 
signals received are ordered in decreasing order of power. For each user, the squared 
moduli of the outputs of the filters adapted to the user connected to the different 
antennae will for example be summed and the sums obtained will be compared. 

In order to simplify the diagram a single stage has been depicted, each stage 

1 5 performing one elimination step. Each stage has at its input a receiver (700k) of the 
first type or of the second type, as illustrated respectively in Figs. 4 and 5. The 
receiver of the first stage receives the L antenna signals and the receiver of a 
subsequent stage receives the L output signals from the previous stage. At the output 
of the receiver, as already described for Fig. 6, a soft or hard decision is made on the 

20 symbols (710 k ) with a demodulation (720k), a deinterleaving and a channel decoding 
(730k), a channel recoding and interleaving (740 k )> a remodulation (750k) and a 
respreading (760k) of the symbols and then a reconstruction filtering (770k)» The 
filters (770 k ) are identical to the filters (670 k ) already described. The L signals output 
from (770 k ) are then subtracted at (780 k ) from the L input signals of the stage k in 

25 order to supply the L input signals of the stage £+1. Thus the contribution of the 
different users are successively eliminated, commencing with those which create the 
greatest interference. The outputs of the reconstruction filters (770 k ) of the different 
stages \..K can be used, like those of the reconstruction filters (670 k ) of Fig. 6, in 
order to regenerate, for each user k, its contribution to the antenna signals, in the form 

30 of an L-tuplet (xejt). The contributions of the different users are subtracted one by 
one by the adders (780 k ) of the different stages. 

The comments made here concerning the initialisation of the coefficients for the 
first iteration apply here for the sequences of the different users. In other words, at 
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each sequence of a user the coefficients are initialised by (18) or (20) if the receiver 
used is of the first type and by (19) or (21) if the receiver used is of the second type. 

In addition, it can be envisaged causing a serial elimination phase (Fig. 7) to be 
followed by a parallel elimination phase (Fig. 6). More generally, the channels can be 
5 grouped together in "packets" corresponding to users of the same spatial region or 
adjacent powers, the interference between packets can be dealt with first of all, and 
then the interference between users of the same packet next. The inter-packet 
treatment can be of a type (serial or parallel) different from the intra-packet treatment. 

10 Fig. 8a depicts a multi-user detection device according to a fourth embodiment 

of the invention. 

This embodiment extends the parallel elimination mode in the sense that it now 
includes a "multi-equalisation" function which will be explained below. The 
reference numbers indicated are the translations of those of Fig. 6 when they 

15 correspond to identical functions. The device depicted in Fig. 8a nevertheless differs 
from the one in Fig. 6 in that it comprises a linear filter (805) upstream of the decision 
modules (81 Ok). The purpose of this filtering is to take into account the outputs of the 
receivers of the first type or second type (800k) at different sampling times in the 
elimination of the interference. The filter (805) is a transverse linear time filter which 

20 fulfils a role equivalent to that of a conventional transverse linear equaliser when a 
single user is considered. Unlike the latter, however, the filter (805) is a matrix filter 
which operates on size vectors K and therefore uses matrices of size KxK instead of 
scalar coefficients, the number of such matrices corresponding to the number of filter 
time slots. The diagonal coefficients of the different matrices allow equalisation 

25 proper whilst the coefficients outside the diagonals make it possible to take into 
account the time responses of the different channels in the rejection of the interference 
between users. Just like its single-user homologue the transverse filter (805) can be 
optimised by means of a known technique such as that of zero forcing (ZF) or 
minimum mean square error (MMSE). 

30 Fig. 8b depicts in an equivalent manner the device of Fig. 8a, showing the vector 

processing character of the different modules (800) to (880), each module respectively 
including the modules (800i v . ; 800 K ) 5 ...(880 1> ...,880 K ) of Fig. 8a (for example, the 
module (800) is a battery of receivers of the first or second type). The matrix filter 
(805) of course remains unchanged. 



19 



Fig. 9 depicts a multi-user detection device according to a fifth embodiment of 
the invention. 

This embodiment transposes the technique of equalisation by decision feedback 
5 (DF) in the multi-user context. The modules (900), (910),... (980) are modules 
identical to the modules (800), (8 10),... ,(880) of Fig. 8, each module operating, as 
stated, on a vector with K components. The device comprises a precursor transverse 
filter (905) identical to the filter (805) and a postcursor filter (907) filtering the vector 
of the symbols $\ at the output of the module (950). The output of the postcursor 

10 filter is subtracted (906) at the input of the decision module (910). The assembly 
consisting of the transverse filter (905) and the feedback filter (907) can, there too, be 
optimised according to a known technique, ZF or MMSE for example. 

The postcursor filter (907) must be causal so as to ensure that there is no need to 
subtract the contribution of time samples not yet available at the time of subtraction of 

15 the current sample. This may notably impose constraints on the length of the 
interleaving. Alternatively, the postcursor filter (907) will filter the symbol vector $k 
at the output of the decision module (910) which will make it possible to relax the 
causality constraint, at the cost however of abandoning the benefit of channel 
decoding in the evaluation of the inter-symbol interference by the postcursor filter. 

20 The embodiments illustrated in Figs. 8a, b and 9 both use a "multi-equalisation" 

derived from single-user equalisation symbol by symbol. It is however clear to a 
person skilled in the art that the invention can also be extended to the family of 
techniques of equalisation with estimation of symbol sequences, namely MLSE 
(Maximum Likelihood Sequence Estimation), DDFSE (Delayed Decision Feedback 

25 Sequence Estimation), RSSE (Reduced State Sequence Estimation) etc. All these 
techniques can also be transposed by means of the invention to the multi-user 
situation. 

Although the invention has been essentially illustrated, for reasons of 
convenience, by means of modules executing certain functions, it is clear however that 
30 all or only some of these functions can be executed by a processor programmed for 
this purpose. 



